Estimate of the charmed hybrid meson spectrum from quenched lattice QCD 
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We compute from quenched lattice QCD the ground state masses of the charmed hybrid mesons 
ccg, with exotic quantum numbers J PC = 1 _ + , +_ and . The hybrid meson spectrum has 
never been provided by lattice simulations due to the difficulties to extract high gluonic excitations 
from noise. We employ improved gauge and fermion actions on the anisotropic lattice, which 
reduce greatly the lattice artifacts, and lead to very good signals. The data are extrapolated to 
the continuum limit, with finite size effects under well control. For 1 ^ and +_ hybrid mesons, 
the ground state masses are 4.405(38) GeV and 4.7f4(52) GeV. We predict for the first time from 
lattice QCD, the ground state mass of __ to be 5.883(146) GeV. 

PACS numbers: 12.38.Gc, 12.39.Mk 



A hybrid (exotic) meson qqg is a bound state of quark 
q, anti-quark q and excited gluon g. The excited gluon 
makes quantum number of the bound state to be 1 ^, 

H or 0~~, , inaccessible to qq mesons in the quark 

model. The existence of hybrid mesons is one of the 
most important predictions of quantum chromodynamics 
(QCD). 

So far no signal for heavy exotic hybrid mesons 
has been experimentally observed, though a number 
of potential candidates for light hybrid mesons were 
suggested 0, 0- Fortunately, this situation may change, 
due to rapid development of new experiments, for ex- 
ample PEP-II(Babar), KEKB (Belle), 12 GeV Jefferson 
LabU EL upgraded CLEO-c detector @, and new BES3 
detector 6J. Especially, 12 GeV Jefferson Lab, CLEO-c 
and BES3 will present well needed and more reliable data 
for the charmonium spectrum, including hybrid mesons. 

The most reliable technique for computing hadron 
spectroscopy is lattice gauge theory. It is a non- 
perturbative approach based on first principle QCD. Of 
course, the lattice approach is not free of systematic er- 
rors. The discretization errors in the Wilson gluon and 
quark actions are the most serious ones. These errors 
are smaller only at very small bare coupling, and very 
large lattice volume is required to get rid of finite size ef- 
fects. The idea of Symanzik improvement |7j is to add new 
terms to the Wilson actions to reduce the lattice spacing 
errors. In combination with tadpole improvement [8|, the 
Symanzik program has recently led to great success in 
approaching the continuum physics on very coarse and 
small lattices. Simulations on anisotropic lattices help 
getting very good signal in spectrum computations. 

There have been many quenched lattice calculations [9|, 
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EimQ[I!ElEEEE[I3 of the 1-+ or 0+- hy- 
brid meson masses, in either light quark or heavy quark 
sector. The mass estimates for the light hybrid mesons 
might still have some uncertainties, because those sim- 
ulations are still far from the chiral regime. The inclu- 
sion of dynamical quarks is still very preliminary^^], due 
to limited computing resources. A recent review can be 
found in Ref.|20|. It has been a long standing puzzle for 
the 0~~ hybrid mesons |Tl|j : no clear signal has ever been 
found, which might be due to the fact that the gluon is 
highly excited. 

As for heavy quarks, special considerations have 
to be taken. Currently, non-relativistic lattice QCD 
(NRQCD), and relativistic heavy quark (Fermilab), and 
anisotropic relativistic approaches are the leading meth- 
ods. Let a denote the lattice spacing. The NRQCD 
method plll^ is applicable for am q > 1 and works well 
for very heavy quarks, especially for the spin-independent 
bb system; however the continuum limit is problem- 
atic because of the condition am q > 1; it is difficult 
to include relativistic corrections and radiative correc- 
tions, leading to breaking down of this method for the 
cc svstem|23j. The relativistic (Fermilab) approach to 
quarks 24J works for both light quarks and heavy quarks; 
Up to 0(a 2 ), the fermionic action is equivalent to the 
standard Sheikholeslami-Wohlert(SW) action|2^| on an 
isotropic lattice; however, to get rid of the 0(a) error 
all coefficients in the fermionic action are required to be 
mass-dependent. The anisotropic relativistic approach 
to quarks 26, 27], which is used in this paper, generalizes 
the Fermilab approach to anisotropic lattice. This im- 
proved quark action has been successfully app lied to the 
computation of the charmonium spectrum|27l l28| . which 
agree very well with experiments. 

To investigate gluonic excitations in hadrons, addi- 
tional improvement of the gluon action would certainly 
help getting better signals. The first attempt was made 
in Ref. ^t|j where the ground state masses of 1 h hybrid 
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mesons in the light quark and charm quark sectors were 
computed, by combining the improved gluon act ion [2^ 
and a simplified relativistic fermionic actional on the 
anisotropic lattice. However, the quark masses were far 
away either from the chiral limit or from the charm quark 
regime. The statistics were low, and finite size effects and 
lattice spacing errors were not analyzed. 

In this letter, we estimate the ground state masses of 
1 ^, -1 , and exotic mesons in the charm quark 
sector, employing lattice QCD with tadpole improved 
gluon|2jJ and quark 20. 27] actions on the anisotropic 
lattice. We get significantly improved signals for these 
particles, in particular for the particle for the first 
time. 

Our simulation parameters are listed in Tab. D At 
each P = 6/g 2 , three hundred independent configurations 
were generated with the improved gluonic action |2{j. 
Two hundred configurations are the minimum for obtain- 
ing stable results. We input two values of bare quark 
mass m q o and then compute quark propagators using 
the improved quark action |27^. and the hybrid meson 
correlation function using the operators 1 = p (£> B, 



o+- = 

in Ref. 



m 



P) ® B, 0p~ = o x (P) ® E and Ol 
as they give the best signals. 
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FIG. 2: Extrapolation of the charmed 1 + , + , and 
hybrid meson masses to the continue limit. 
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FIG. 3: Ratio of splittings Rh = AM (IH -IS)/ AM (I 1 Pi - 
IS) against a 2 s . The straight line is the extrapolation to con- 
tinuum limit. 



FIG. 1: Effective mass of the hybrid meson for /3 = 3.0 
and atrriqo = 0.100. The solid line is the fitted result, ranging 
from ti = 6 to tf = 12 with \ 2 /d.o.f. = 0.4326 and confidence 
level=0.7620. 

Figure ^ shows an example of the effective mass plot 
]n(C(t)/(C(t + 1)) of the 0— hybrid, where C(t) is the 
correlation function between the hybrid operators. We 
obtained the effective mass atm, with the fit range cho- 
sen according to optimal confidence level and reasonable 
X 2 /d.o.f. 

We then interpolated the data to the charm quark 
regime using (m„(K t — * K charm^ _|_ 3 TOp ( Ki 
/ff orm ))/4 - M(lS) exp (m(^ c ) e;cp + 

3m(J/tfj) exp )/A = 3067. 6MeV, where the right hand side 
is the experimental value for the 15 charmonium. The 



results for the charmed hybrid meson masses are listed 
in Tabs. |H]for the ground state. It is also important to 
check whether these lattice volumes are large enough. 
We also did simulations on 8 3 x 48 and 12 3 x 48 at 
/3 = 2.6, 12 3 x 36 at /3 = 2.8, and 16 3 x 48 at (3 = 3.0, 
but here we just list the results from the largest volume. 
When the spatial extent is greater than 2.2fm, the finite 
volume effect on the _ ~ mass is less than 0.1% for the 
ground state. 

The l 1 P\ — IS charmonium splitting was chosen to de- 
termine the lattice spacing, because it is roughly indepen- 
dent of quark mass for charm and bottom sectors, and 
the experimental value AM(l 1 Pi - lS) exp = 457MeV 
was well measured. Here we used \ Cl meson (1 ++ ) mass 
for the P-wave and m(7y c )/4 + 3m(J/i/ , )/4 for the S-wave. 
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TABLE I: Simulation parameters at largest volume. We employed the method in Ref.|27| to tune these parameters, Kt and 
Us for the quark action. The last two columns are about the spatial lattice spacing and the lattice extent in physical units, 
determined from the IP — IS charmonium splitting. 
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TABLE II: Charmed hybrid meson spectrum for the ground state. The results in the continuum limit (/3 = oo) were obtained 
by: (i) directly extrapolating the data to a 2 s — ► 0; and (ii) using the ratio of splittings Rh, as described in the text. 



The results for the spacial lattice spacing a s at different 
(3 are listed in Tab. [I] They are consistent with those 
from the heavy quark potential|29l|. 

Figure [3 shows the charmed 1 h , H , hybrid 
meson masses as a function of a 2 s . Other hybrids have 
a similar behavior, indicating the linear dependence of 
the mass on a 2 .. The spectrum in the continuum limit is 
obtained by linearly extrapolating the data to a 2 —* 0, 
as also listed in Tabs. [H] 

We also computed the splitting of 1H — IS and the 
ratio R H = AM(1H - IS)/AM{1 1 P 1 - IS), where 1H 
stands for the ground state of a charmed hybrid meson. 
Its dependence on a 2 and extrapolation to the continuum 
limit are shown in Fig. |3| The last line of Tab. [n] also 
lists the hybrid meson masses for the ground state, using 
the following equation: 

lim M(IH) = M(lS) exp 

a s~ *0 

+ AM(1 1 P 1 - lS) exp x lim R H . 

This method was claimed to be better because the 
splitting between a hybrid and the IS" state is rather 
insensitive to the imperfect tuning of AM(l 1 Pi — IS) 
and M(1S). However, as seen in Tab. [H] the results 
from two different methods agree very well. 

One source of systematic errors in our calculation is 
the quenched approximation. Although full QCD simu- 
lations will remove this unknown error, quenched approx- 



imation in some areas[2ll including the hybrids j^J, 
continues to play an important role. The findings in Refs. 

0] indicate that the effects of dynamical quarks 
on light hybrids and bottomed hybrids are very small. To 
have full relevance of the charmed hybrids to experiment, 
simulations with dynamical quarks, although extremely 
expensive to achieve high statistics, might be helpful to 
see whether the quenching error is under control. Nev- 
ertheless, our results are a very important step for com- 
parison with future dynamical simulations. 

Simulations on the anisotropic lattice with both gluon 
action and improved quark action improved lead to the 
first observation of the clear signal for the hybrids. 
We believe that our findings are useful to experimental 
search for these new particles, predicted by QCD. 
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